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ABSTRACT 

Studying the nature and origin of the intergalactic magnetic field (IGMF) is 
an outstanding problem of cosmology. Measuring Faraday rotation would be a 
promising method to explore the IGMF in the large-scale structure (LSS) of the 
universe. We investigated the Faraday rotation measure (RM) due to the IGMF 
in filaments of galaxies using simulations for cosmological structure formation. 
We employed a model IGMF based on turbulence dynamo in the LSS of the 
universe; it has an average strength of (B) ~ 10 nG and a coherence length of 
several x 100 h~ x kpc in filaments. With the coherence length smaller than path 
length, the inducement of RM would be a random walk process, and we found 
that the resultant RM is dominantly contributed by the density peak along line 
of sight. The rms of RM through filaments at the present universe was predicted 
to be ~ 1 rad m~ 2 . In addition, we predicted that the probability distribution 
function of |RM| through filaments follows the log- normal distribution, and the 
power spectrum of RM in the local universe peaks at a scale of ~ 1 h~~ l Mpc. 
Our prediction of RM could be tested with future instruments. 

Subject headings: intergalactic medium — large-scale structure of universe - 
magnetic fields — polarization 
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Introduction 



Th e intergalactic m edium (IGM) contains gas, which was heated mostly by cosmological 
shocks (IRyu et al.ll2003l ). along with dark matter; the hot gas with T > 10 7 K is found inside 
and around clusters/groups of galaxies and the warm- hot intergalactic medium (WHIM) 
with T = 10 5 — 10 7 K resides mostly in filaments of gal axies, while lower temperature gas 



is dis tributed mostly as sheetlike structures or in voids (ICen fc Ostrikerl Il999l ; iKang et al. 



20051 ). As the gas in the interstellar medium, the gas in the intracluster medium (ICM) and 
filaments is expected to be permeated with magnetic fields. Measuring Faraday rotation, 
the rotation of the plane of linearly-polarized light due to the birefringence of magneto-ionic 
medium, has been one of a few methods to explore the intergalactic magnetic field (IGMF). 

Observational exploration of the IGMF using Faraday rotation me asure (RM) was 
start ed with the investigation of the intracluster magnetic field (ICMF) (see ICarilli &: Taylor 
20021 . for a review). An RM study of the Coma cluster, for instance, reve aled the ICMF o f 
the strength of order ~ fiG for the coherent length of order ~ 10 kpc (IKim et al.l Il990l ). 
For Abell clusters, the RM of typically ~ 100 — 200 rad m~ 2 was obseryed, in dicating an 
average strength of the ICMF to be ~ 5-10 /iG rtciarke et allboOll ; Idarkel 120041 ) . RM maps 



of clusters were analyzed to study the power spectrum of turbulent magnetic fields in the 
ICM; for instance, a Kolmogorov-like spectrum with a bending at a few kpc scale was found 
in the cooled core region of the Hydra cluster (IVogt fc En61mll2005l ). and spectra consistent 
with the Kolmogoroy spectrum were reported i n the wider ICM for t he Abell 2382 cluster 
( IGuidetti et al.ll2008l ) and for the Coma cluster (IBonafede et al.ll2010l ). 



The nature of the IGMF in filaments, on the co ntrary, remains largely unknown, because 
the study of RM outside clusters is still scarce (e.g-. Ku et al.ll2006l ); detecting the RM due to 



the IGMF in filaments is difficult with current facilities, and also removing the galactic fore- 
ground is not a trivial task. The next generation radio interferometers including the Square 
Kilometer Array (SKA), and upcoming SKA pathfinders, the Australian SKA Pathfinder 
(ASKAP) and the South African Karoo Array Telescope (MeerKAT), as well the Low Fre- 
quency Array (LOFAR), however, are expected to be used to study the RM. Particularly, 
the SKA could measure RM for ~ 10 8 polarized extragalactic sources acr oss the sky with an 
average spacing of ~ 60 arcsec between lines of sight (LOS's) (see, e.g., ICarilli fc Rawlings 
2004 ; iKrause et al.ll2009l . and references therein), enabling us to investigate the IGMF in the 
large-scale structure (LSS) of the universe. 

Att empts to t heore ticall y predict the RM due to the IGMF have been made: for in- 



stance, IRyu et al.l (119981 ) and iDolag et al.l (120051 ) used hydrodynamic sim ulations for cosmo- 
l ogica l structure formation to study RM in the LSS, and more recently iDubios fc Teyssier 
(120081 ) used MHD simulations to study RM for clusters. However, the properties of the 
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IGMF, especially in filaments, such as the strength and coherence length as well as the spa- 
tial distribution, are largely unknown, hindering the theoretical study of RM in the LSS of 
the universe. 



Recently, iRyu et al.l ( 120081 ) proposed a physically motivated model for the IGMF, in 
which a part of the gravitational energy released during structure formation is transferred to 
the magnetic field energy as a result of the turbulent dynamo amplification of weak seed fields 
in the LSS of the universe. In the model, the IGMF follows largely the matter distribution 
i n the cosmic web and the strength is predicted to be (B) ~ 10 nG in filaments. ICho fc Ryu 
( 120091 ) studied various characteristic length scales of magnetic fields in turbulence with very 
weak or zero mean magnetic field, and showed that the coherence length defined for RM is 
3/4 times the integral scale in the incompressible limit. They predicted that in filaments. 



the co herence length for RM would be a few x 100 h 1 kpc with the IGMF of IRyu et al. 
( 120081 ) and the RM due to the magnetic field would be of order ~ 1 rad mT 2 . 



In this paper, we study RM in the LSS of the universe, focusing on RM through fila- 
m ents, using simul ation s for cosmological structure formation along with the model IGMF 
of IRyu et al.l (120081 ) and lCho fc Ryul (120091 ) . Specifically, we present the spatial distribution, 
probability distribution function (PDF) and power spectrum of the RM, and discuss the 
prospect of possible observations of the RM. In sections 2 and 3, we describe our model and 
the results. Discussion is in Section 4, and Summary and Conclusion follows in Section 4. 



Model 



To investigate RM in the LSS of the universe, we used structure formation simulations 
for a concordance ACDM universe with the following values of cosmological parameters: 
Qbm = 0.043, Q DM = 0-227, Q\ = 0.73, h = #o/(100 km s" 1 Mpc" 1 ) = 0.7, n = 1, and 



erg = 0.8 (same as in lRyu et al.ll2008l). Thev were pe rformed using a particle-mesh/Eulerian, 
cosmological hydrodynamic code (IRyu et al.l 119931 ). A cubic region of comoving volume 
(100 /i _1 Mpc) 3 was reproduced with 512 3 uniform grid zones for gas and gravity and 256 3 
particles for dark matter, so the spatial resolution is 195 h^ 1 kpc. Sixteen simulations with 
different realizations of initial condition were used to compensate cosmic variance. 



For the IGMF, we employed the model of IRyu et al.l ( 120081 ); it proposes that turbulent- 
flow motions are induced via the cascade of the vorticity generated at cosmological shocks 
during the formation of the LSS of the universe, and the IGMF is produced as a consequence 
of the amplification of weak seed fields of any origin by the turbulence. Then, the energy 
density (or the strength) of the IGMF can be estimated with the eddy turnover number and 
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the turbulent energy density as follow: 



£b 



^eddy 



^turb- 



Here, the eddy turnover time is defined as the reciprocal of the vorticity at driving scales, 
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eddy 
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(u = V x v), and <fi is the conversion factor from turbulent to magnetic 
energy that depends on the eddy turnover number t/t eddy . The eddy turnover number was 
estimated as the age of universe times the magnitude of the local vorticity, that is, t age oj. The 
local vorticity and turbulent energy density were calculated from simulations for cosmological 
structure formation described above. A functional form for the conversion factor was derived 
from a separate, incompressible, magnetohydrodynamic (MHD) simulation of turbulence 
dynamo. For the direction of the IGMF, we used that of the passive fields from simulations for 
cosmological structure formation, in which weak seed fields were evolved passivel y, ignoring 
the back-reaction, along with flow motions (IKulsrud et al.l 1 19971 ; iRyu et al.lll998[ ). 



In our model, a s seed magnetic fi elds, we took the ones generated through the Biermann 
battery mechanism (jBiermann!ll950l ) at cosmological shocks. There are, on the other hand, 
a number of mechanisms that have been suggested to create seed fields in the early universe. 
Besides various inflationary and string theory mechanisms, the followings include a partial 
list of astrophysical mechanisms. At co smological shocks, in a d dition, Weibel instability ca n 
operate and produce magnetic fields (IMedvedev et al.l 120061 ; ISchlickeiser k, Shuklal 120031 ) . 
and streaming cosmic rays accelerate d by the s hocks can amplify weak upstream magnetic 
fields via non-resonant growing mode (jBelll 120041 ) . In addition, for ins tance, galactic outflows 
during the starburst phase of galactic evolution ( Donnert et al. 2009h and the return c urrent 
induced by cosmic-rays produced by Supernovae of first stars (jMiniati fc Belli 120101 ) were 
suggested to deposit seed fields. We point, however, that in our model the IGMF resulting 
from turbulent amplification should be insensitive to the origin of seed fields. 

The spa ti al dis tribution of the strength of the r esulting IGMF is shown in Figure 4 
of IRyu et al.l (120081 ) and Figure 1 of IRyu et al.l (120101 ) . It is very well correlated with the 



distribution of matter. The average strength of our model IGMF for the WHIM with 10 5 < 



T < 10 7 K in filaments is (B) ~ 10 nG, {B 2 ) 1/2 
((pfi) 2 ) 1 / 2 /(p 2 ) V2 ~ a few x 0.1 nG. 



a few x 10 nG, (pB)/(p) ~ 0.1 fiG, or 



3. Results 



We calculated RM, defined as Ax/AA 2 (% is the rotation angle of linearly-polarized 
light at wavelength A), in the local universe with z = along a path length of L = 100 hT 1 
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Mpc, which is the box size of structure formation simulations. Figure 1 shows the resulting 
RM map of (28 /i" 1 Mpc) 2 area in logarithmic and linear scales. RM traces the large-scale 
distribution of matter, and we see two clusters and a filamentary structure containing several 
groups. Through the clusters, groups, and filament in the field, RM is roughly ~ 100, ~ 10, 
and ~ 1, respectively, while RM through sheets and voids is much less. The bottom panel 
of Figure 1 shows the mixture of positive and negative RM, reflecting the randomness of 
magnetic fields in the LSS. 

With the coherence length of magnetic fields for RM (see Discussion) expected to be 
smaller than the path length which should be a cosmological scale, the inducement of RM is 
expected be a random walk process. Figure 2 shows the distributions of RM as well as other 
quantities along a few LOS's through filaments; it confirms that the inducement of RM is 
indeed a random walk process. However, we note that the resulting RM is dominated by the 
contribution from the density peak along LOS's. 

To quantify RM in the LSS of the universe, we calculated the probability distribution 
function (PDF) of |RM| for 512 2 x 3 x 16 (projected grid zones x directions x runs) LOS's. 
Figure 3 shows the resulting PDF through the LOS's of different ranges of the mean temper- 
ature weighted with X-ray emissivity, T X - The figure also shows the fitting to the log- normal 
distribution, 



PDF(log 10 |RM|) 



V2 



: exp 



(log 10 |RM|-//)' 
2a 2 



(2) 



finding that the PDF closely follows the log-normal distribution. We also calculated the 
root mean square (rms) of RM, (RM) rms ; note that the mean of RM, (RM), is zero for our 
IGMF. Through the WHIM, which mostly composes filaments, (RM) rms = 1.41 rad m~ 2 . 
This agrees well with th e value predicted w ith the mean strength and coherence length of 
the IGMF in filaments by lCho &: Ryul (120091 ) . However, this is an order of magnitude smaller 
than the values of |RM| toward the Hercules and Perseus-Pisces superclusters reported in 
Xu et al.l (120061) . Th e diffe rence is mostly due to the mass-weighted path length; the value 
quoted by Ku et al.l (120061 ) is about two orders of magnitude larger than ours. Through 
the hot gas with 10 7 < T < 10 8 K, (R M) rms = 108 rad m" 2 , whi ch is in good agreement 
with RM observations of galaxy clusters (IClarke et al.ll200ll ; IClarkd 12004 ). Through the hot 
gas, however, we found RM of up to > 1000 rad m~ 2 . This should be an artifact of limited 
resolution (see Discussion). So the values for the hot gas in our work should not be taken 
seriously. 

Finally, we calculated the two-dimensional power spectrum of RM on 3 x 16 (directions x 
runs) projected planes; Prm(^) ~ |RM(/c)| 2 fc, where RM(fc) is the Fourier transform of 
RM(f) on planes. Figure 4 shows the resulting power spectrum along with the power spectra 
of electron density, magnetic fields, and the curl component of flow motions, v cur b which 
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satisfies the relation V x v CUT \ = V x v . The power spectrum of RM peaks at k ~ 100, 
which corresponds to ~ 1 h^ 1 Mpc. Cosmic variance is not significant around the peak, 
although it is larger at smaller k, as expected. The power spectrum of RM reflects the 
spatial distributions of electron density, n e , and LOS magnetic field, B\\. The power spectra of 
projected n e and projected B\\, have peaks at ~ 3 h^ 1 Mpc and ~ 1.5 h^ 1 Mpc, respectively. 
The shape of the power spectrum of RM follows that of projected B\\ rather than that of 
projected n e , implying that the statistics of RM would primarily carry the statistics of the 
IGMF. 



4. Discussion 



Our results depend of RM on the strength and coherence length of the IGMF. We 
employed a model where the strength of the local IGMF was estimated based on turbulence 
dynamo, while the direction was gripped from structure formation simulations with passive 
fields (see Section 2). In principle, if we had performed full MHD simulations, we could have 
followed the amplification of the IGMF through turbulence dynamo along with its direction. 
In practice, however, the currently available computational resources do not allow a numerical 
resolution high enough to reproduce the full development of MHD turbulence. Since the 
numerical resistivity is larger than the physical resistivity by many orders of magnitude, the 
growth of magnetic fields is expected to be saturated before dynamo action becomes fully 
operative (see, e.g., Kulsrud et al.lll997| ). In such situation, the state of magnetic fields in 
full MHD, including, for instance, the power spectrum, is expected t o mim ic that of passive 
fields. This is the reason why we adopted the model of iRyu et al.l (120081 ) to estimate the 
strength of the IGMF, but we still used passive fields from structure formation simulations 
to model the field direction. 



The validity of our model IGMF was checked as follows: 

1) In MHD turbulence, the distribution of magnetic fields, including the direction, is ex- 
pected to correlate with that of vorticity, since magnetic fields and vorticity are described by 
similar equations except the baroclinity term in t he equation for vorticity (if dissipative pro- 
cesses are ignored) (see, e.g.. Kulsrud et al.lll997| ). Such a correlation can be clearly seen in 
Figure 5, in which we depicts the distributions of our IGMF and vorticity in two-dimensional 
slices. 

2) Full MHD turbulence simulations suggest that the peak of magnetic field spectrum occurs 
~ 1/2 of the energy injection scale, or the peak scale of velocity power spectrum, at satu- 
ration : in the linear gro wth stage, the peak scale of magnetic field spectrum grows as ~ t 1 ' 5 



or so 



icho fc Rvull2009l l With our model IGMF, the peak scale of magnetic field spectrum 
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is ~ 1 h^ 1 Mpc (the third panel of Figure 4); on the other hand, the curl component of 
flow motions has the peak of power spectrum at ~ 4 h~ l Mpc (the bottom panel of Figure 
4). That is, the peak scale of magnetic field spectrum is ~ 1/4 of the energy injection scale 
in our model IGMF. By considering the t urbulence in the LSS of the universe has not yet 



reached the fully saturated stage (see, e.g-. lRyu et al.ll2008l ). the ratio of the two scales seems 
to be feasible. 

These suggest that our model IGMF would produce reasonable results, although eventu- 
ally it needs to be replaced with that from full MHD simulations for cosmological structure 
formation when computational resources allow such simulations in future. 

Apart from our model for the IGMF, the finite numerical resolution of simulations 
could affect our results. The average strength of our model IGMF is (B) ~ a few iiG in 



cluste rs/groups, ~ 0.1/iG around clusters/groups, and ~ 10 nG in filaments. iRyu et al. 



(120081 ) tested the numerical convergence of the estimation. With simulations of different 
numerical resolutions for cosmological structure formation, it was shown that (B) of our 
model IGMF for the WHIM with 10 5 < T < 10 7 K would approach the converg ence value 



within a factor ~ 2 — 3 at the resolution of 512 3 grids (see Figure S5 of SOM of IRyu et al. 



(12008( 1 ). 



It is, on the other hand, rather tricky to assess the effect of finite resolution on the 
coherence length of our model IGMF, because the definition of coherence length for RM is 
not completely clear and the estimation of coherence length, for instance, for the filament 
IGMF alone is not trivial. We tried to quantify coherence length in the following three ways: 
1) We directly calculated the coherence length of B\\, that is, the length with the same sign 
of Bn, along LOS's. Figure 6 shows the PDF of the resulting coherence length through the 
WHIM, which composes mostly filaments. It peaks at the length of 3 zones corresponding 
to 586 hr x kpc. 2) We calculated 3/4 times the integral scale, 

3 jP* D {k)/k dk 

4 X 2n jPl D {k)dk ' (3) 

which is the coherence length defined for RM in the incompressible limit (see Introduction), 
for the IGMF inside the whole computational box of (100 ^i _1 Mpc) 3 volume. Here, P^ D (k) 
is the three-dimensional power spectrum of magnetic fields (the third panel of Figure 4). 
We found the value to be ~ 800 h^ 1 kpc for our model IGMF. 3) We also calculated the 
largest energy containing scale in the whole computational box, which is the peak scale of 
kPl D (k) (not shown). It is ~ 900 h 1 kpc for our model IGMF. Note that the latter two 
values include contributions from the IGMF in filaments as well as in clusters, sheets, and 
voids. All the three scales are comparable. These length scales are ~ 3 to 5 times larger 
than the grid resolution of our simulations, 195 h^ 1 kpc. 
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Cho Sz Ryul (120091 ) studied characteristic lengths in incompressible simulations of MHD 
turbulence (see Introduction); based on it, they predicted that the coherence length for RM 
would be a few x 100 h^ 1 kpc in filaments, while a few x 10 h^ 1 kpc in clusters. With our 
grid resolution of 195 h~ l kpc, the coherence length of the IGMF in clusters should not be 
resolved and so our estimation of RM for clusters should be resolution-affected, as pointed 
in Section 3. On the other hand, while the predicted coherence length for the IGMF in 
filaments is still larger than the grid resolution, the estimated coherence length of B» for the 
WHIM is a couple of times larger than the prediction for filaments. It could be partly due to 
the limited resolution in our simulations. However, as noted in Section 3, RM is dominantly 
contributed by the density peak along LOS's (Figure 2). 

The above statements indicate that our estimate of the RM through filaments is expected 
to have uncertainties, especially due to limited resolution of our simulations; the error in our 
estimation could be up to a factor of several. 



5. Summary and Conclusion 



We studied RM in the LSS of the universe, focusing on RM through filam ents; simu- 



l ation s for cosmological struc ture formation were used and the model IGMF of iRyu et al. 



(120081 ) and I Cho &: Ryul (120091 ) based on turbulence dynamo was employed. Our findings are 
summarized as follows. 1) With our model IGMF, the rms of RM through filaments at the 
present universe is ~ 1 rad m -2 . 2) The PDF of |RM| through filaments follows the log- 
normal distribution. 3) The power spectrum of RM due to the IGMF in the local universe 
peaks at a scale of ~ 1 h^ 1 Mpc. 4) Within the frame of our mode IGMF, we expect that 
the uncertainty in our estimation for the rms of RM through filaments, due to the finite 
numerical resolution of simulations, would be a factor of a few. 

We note that our model does not include other possible contribu tions to the IGMF, fo r 
instance, that from galactic black holes (AGN feedbacks) (see, e.g., Kronberg et al.ll200ll ). 
So our model may be regarded as a minimal model, providing a baseline for the IGMF. 
With such contributions, the real IGMF might be somewhat stronger, resulting in somewhat 
larger RM. 

It has been suggested that future radio observatories such as LOFAR, ASKAP, Mee rKAT 



and SKA could detect the extragalactic RM of ~ 1 rad m~ 2 we predict (see, e.g., [Beck 
20091 ). However, it is known that the typical galactic foreground of RM is a few tens 
and of order ten rad m~ 2 in the low and high galactic latitudes, respectively (see, e.g., 



Simard-Normandin &: Kronberg I Il980l ). So the detection of the extragalactic RM of 
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1 rad m -2 or so could be possible only after the galactic foreground is removed. We note, 
however, that with filaments at cosmological distance, the peak of the power spectrum of 
the RM due to the IGMF in filaments would occur at small angular scales; for instance, for 
a filament at a distance of 100 h^ 1 Mpc, the peak would occur at ~ 0.5 degree or so. This is 
much smaller than the expected angular scale of the peak of the galactic foreground, which 



would be around tens degree (see, e.g., iFrick et al.l 1200 ll ). Then, it would be plausible to 



extract the signature of the RM of ~ 1 rad m~ 2 due to the IGMF in filaments. We leave 
this issue and connecting our theoretical prediction of RM in the LSS of the universe to 
observation for future studies. 

TA was supported in part by National Research Foundation of Korea (R0 1-2007-000- 
20196-0). DR was supported in part by National Research Foundation of Korea (K20901001400- 
09B1300-03210) 
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Fig. 1. — RM map of (28 /i _1 Mpc) 2 area in the local universe of depth (path length) of 
L = 100 hT 1 . Top and bottom panels show the map in logarithmic and linear scales, 
respectively. 
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RM=2.91 RM= 0.886 RM= 0.316 
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Fig. 2. — Profiles of local RM, LOS magnetic field, B», electron density, n e , and gas tem- 
perature, T, along a few LOS's through filaments. The arrows indicate the sign of local RM 
and B\\. 
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Fig. 3. — PDF of |RM| through LOS's of different ranges of the mean temperature weighted 
with X-ray emissivity, Tx- Thin lines, thick lines (red or black), and thick lines (blue or 
gray) show the PDFs from 16 independent runs, their average, and the best-fit to the log- 
normal distribution, respectively. The values of fitting parameters and the rms of RM are 
also shown. 
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Fig. 4. — Two-dimensional power spectra of RM, projected electron density (n eiPro j), pro- 
jected line-of-sight IGMF strength (-B|| )Pxo j), and projected curl component of flow motions 
(v cm \) from top to bottom panels, respectively. Thin solid lines show the power spectra for 
3x16 two-dimensional maps, and thick solid lines show their average. The three-dimensional 
power spectra of n e , B, and v cm \, are also shown with thick dashed lines. 
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Fig. 5. — Distributions of the IGMF (left) and vorticity (right) in two-dimensional slices. 
The length of arrows for the IGMF corresponds to x of 10 x_12 G. The color shows the gas 
temperature. 
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Fig. 6. — PDF of the coherence length of B\\ along LOS's (the length with the same sign of 
B\\) through the WHIM. Thin and thick lines show the PDFs for 16 independent runs and 
their average, respectively. 



